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Abstract: 

The existing calculations of the form factors describing the decay B —> peu from QCD 
sum rules have yielded conflicting results at small values of the invariant mass squared 
of the lepton pair. We demonstrate that the disagreement originates from the failure of 
the short-distance expansion to describe the p meson distribution amplitude in the region 
where almost the whole momentum is carried by one of the constituents. This limits 
the applicability of QCD sum rules based on the short-distance expansion of a three- 
point correlation function to heavy-to-light transitions and calls for an expansion around 
the light-cone, as realized in the light-cone sum rule approach. We derive and update 
light-cone sum rules for all the semileptonic form factors, using recent results on the p 
meson distribution amplitudes. The results are presented in detail together with a careful 
analysis of the uncertainties, including estimates of higher-twist effects, and compared to 
lattice calculations and recent CLEO measurements. We also derive a set of “improved” 
three-point sum rules, in which some of the problems of the short-distance expansion are 
avoided and whose results agree to good accuracy with those from light-cone sum rules. 
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1 Introduction 


The interest in the study of semileptonic B decays is mainly due to their importance in 
determining the CKM matrix elements \Vcb\ and \Vub\- Whereas the theoretical analysis 
of both exclusive and inclusive b ^ c transitions is decisively simplihed by an expansion 
in the inverse heavy quark mass, this method is of only little use in b ^ u transitions. 
This is essentially due to the fact that in inclusive channels experimental observation 
is possible only in a small region of phase-space beyond the kinematical threshold for 
charm production, in which the hadron multiplicity is small. Thus, since the theoretical 
description of inclusive decays is essentially based on a parton-model picture, it is not 
very predictive in the experimentally accessible range, cf. ]l[|. It is therefore rather the 
exclusive decay channels, in particular B niu and B —>■ piu, that seem to be more 
suitable for obtaining reliable information on \Vub\- The CLEO collaboration has recently 
presented first experimental results of these branching ratios p|, which, however, are model- 
dependent. 

The decay B —> vrf'z/ has been tackled by several authors using a number of different 
approaches, in particular QCD sum rules H and simulations on the lattice [ 0 , |, |, 

pi]| ; the results are in reasonable agreement. The situation is, however, not that favourable 
in the B ^ p channel. Although here the same methods were applied, the resulting 
predictions for the decay rates are quite different. To illustrate the origin of the problem, 
let us hrst introduce the relevant observables: the hadronic matrix element determining 
the 5 —> p weak transition is 

{p,\\{V - A)^\B) = -i{mB + 

^ tub -|- rrip 


iAAt) , X 2V{t) 

+- - Pp)^ +-^ 


niB + m 


tub + m 




( 1 . 1 ) 


where the four form factors Ai^ 2,3 and V depend on the momentum transfer t = {pB—PpY to 
the leptons; in the limit of vanishing lepton mass A 3 does not contribute to the semileptonic 
decay rate and will not be considered in this paper. A is the polarization of the p, t 
takes values between 0 and 20.3 GeV^. It is precisely this large range of relevant t that 
renders the theoretical description of the form factors in ( p..l|) so difficult. Most quark 


model calculations rely essentially on the pole-dominance picture |]^ or on a nonrelativistic 
description, which yields an exponential increase of the form factors with t [^], which was 
softened in an updated version of the model. Ref. [IR . Lattice calculations are up to now 


limited to small momenta of the hnal state p Q and/or require extrapolations in the 
heavy quark mass. For B —>■ n£u, the possibility to restrict the functional t dependence 
of the single relevant form factor from unitarity with the supplementary input of available 
lattice data at large t was investigated in |TO, but this method is presently not very 


predictive in the p channel, see 

To date, only lattice simulations and QCD sum rules seem to be apt to predict the 
t dependence in nearly the whole physical region. QCD sum rules provide a consistent 
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Figure 1: “Hard” (a,b) and “soft” (c) contributions to the decay form factor, (d): modelling 
of the soft contribution in the QCD sum rule approach. See explanation in the text. 


description of semileptonic D decays |]^ and of the i? —> tt semileptonic form factor 
i, I, I, i- However, there exist conflicting predictions from different types of QCD sum 
rules for B ^ p decays, which differ by a factor 2 in the form factors at maximum recoil 
EM- The aim of this paper is to clarify the origin of this difference and to give updated 
predictions for the form factors, which include in particular recently gained information on 
the structure of p mesons probed at large momentum transfer |T^ . 


At large recoil, the light quark originating from the weak decay carries a large energy 
of order mb/2 and has to transfer it to the soft cloud to recombine to the hnal state 
hadron. The probability of such a recombination depends on the parton content of both 
the B meson and the light meson, the valence conhguration with the minimum number of 
Fock constituents being dominant. The valence quark configuration is characterized by the 
wave function 0(x, k±) depending on the momentum fraction x carried by the quark and on 
its transverse momentum k±. There exist two different mechanisms for the valence quark 
contribution to the transition form factor. The first one is the hard rescattering mechanism 
pictured in Figs. l(a,b), which requires that the recoiling and spectator quarks are at small 
transverse separations. In this case the large momentum is transferred by an exchange of 
a hard gluon with virtuality ~ 0{mb). This contribution is perturbatively calculable in 
terms of the Bethe-Salpeter wave functions at small (~ 1/mb) transverse separations, or 
distribution amplitudes-. 


(j){x) = 


dk‘]_(j){x, k±). 


( 1 . 2 ) 


The second mechanism is the soft contribution, shown schematically in Fig. 1(c). The 
idea is that hard gluon exchange is not necessary, provided one picks up an “end-point” 
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configuration with almost all momentum 1 — x ~ (9(l/mfe) carried by one constituent. 
Since at large scales (p{x) ~ 1 — x [^, the overlap integral is of order dx(j){x) ~ 


0(l/m^). An additional factor comes from the normalization of the heavy initial state, 
so that the hnal scaling law for the soft contribution to the form factor at large recoil is 
1^ Note that the transverse quark-antiquark separation is not constrained 


in this case, which implies that the soft contribution is sensitive to long-distance dynamics. 
To calculate the soft contribution one needs to know the wave function as a function of 
the transverse separation; the simpler distribution amplitude is not enough. 

Hard exclusive processes involving light hadrons receive the same two types of contribu¬ 
tions. There is a major difference, however, in that for light hadrons the soft contribution 
is suppressed by a power of the large momentum (i.e. it is of higher-twist), while for heavy 
meson decays both soft and hard contributions turn out to be of the same power in the 
heavy quark mass B0. PT| . As a result, the soft (end-point) contribution is expected 
to be large and requires quantitative evaluation. 

It was suggested that Sudakov-type perturbative corrections cut off contributions 
of large transverse separations so that the soft contributions might be suppressed. This 
suppression eventually eliminates the soft contribution for very large rri},. At the realistic 
values mb ^ 5 GeV, however, it is unlikely that calculations of this type can provide a 
quantitative description. Indeed, the existing estimates of “hard” contributions typically 
fall short of realistic values of the form factors from model calculations.^ 

Here the QCD sum rules method enters the stage and suggests a nonperturbative 
technique to estimate the necessary convolution integral without explicit knowledge of the 
wave functions. 

There exist two different types of QCD sum rules, which, being similar in spirit, differ 
in the treatment of the light hadron in the hnal state. This is illustrated in Figs. l(c,d). 

The “traditional” sum rules avoid introducing wave functions altogether by considering 
a three-point correlation function with a suitable interpolating current and use dispersion 
relations to extract the contribution of the ground state. The most important nonpertur¬ 
bative effect is then described by the diagram in Fig. 1(d), where the light quark is soft and 
interacts with the nonperturbative QCD vacuum, forming the so-called quark condensate. 
Since quarks in a condensate have zero momentum, it is clear that this diagram yields 
a contribution to the distribution amplitude that is naively proportional to 5(1 — x) and 
remains unsuppressed for mb —> oo. This obviously violates the power counting discussed 
above. The contradiction must be resolved by including the contributions of higher-order 
condensates to the sum rules and subtracting the contribution of excited states. The sup¬ 
pression of the end-point region x —>• 1, which is strictly required by QCD, is thus expected 
to hold as a numerical cancellation between different contributions, which becomes the 
more delicate (and requires more hne-tuning) the more increases. For mb ~ 5 GeV 
a suppression of the quark condensate contribution by a factor IGeV^/m^ ~ 1/25 is 

^ Note that a similar suppression is present for hard exclusive processes involving light hadrons |2^ (in 
which case the soft contribution is additionally suppressed by a power of the large momentum); however, 
there is increasing evidence that soft contributions to, say, the pion form factor remain important at least 
up to « 10 GeV^, see 
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required. This explains why the traditional three-point sum rules become unreliable. 

The light-cone sum rules avoid this problem by rearranging the calculation in such a way 
that nonperturbative effects like the interaction with the quark condensate are included 
in the nonperturbative hadron distribution amplitudes, estimated using additional sum 
rules m- These additional sum rules are written for integrated characteristics of the 
distribution amplitudes like moments, and the correct asymptotic behaviour at the end¬ 
points is included by construction. 

The contribution of a single leading-twist distribution amplitude incorporates an infinite 
set of contributions of condensates of increasing dimension in the standard approach, at 
the cost of retaining pieces with the largest Lorentz spin (lowest twist) only. The remaining 
pieces of condensate contributions are organized in a similar way in the contributions of 
higher-twist distribution amplitudes.^ 

The premium for this rearrangement is that light-cone sum rules have the correct asymp¬ 
totic behaviour in the heavy quark limit, but the snag is that the present knowledge of 
higher-twist distribution amplitudes is incomplete, so that not all known nonperturbative 
corrections of the standard approach can be included. One should expect that these two 
approaches provide complimentary descriptions of B decays, with their own advantages 
and disadvantages. 

Note that the problem with three-point sum rules originates from the constraint of the 
distribution amplitude convolution integral to the end-point region; this makes the answer 
very sensitive to the precise shape of the distribution amplitude rather than its integrated 
characteristics. For the decay form factors at small recoil there is no such restriction, and 
the classical QCD sum rule concept of taking into account nonperturbative effects as the 
contributions of long-wave vacuum fluctuations (condensates) classified by their dimension 
is adequate. Thus, at small recoil, both types of the sum rules should yield similar results 
and the spread of their predictions characterizes the accuracy of the method. At large 
recoil one should rely on the light-cone sum rules. 

It is worthwhile to add that contributions of hard rescattering can be consistently 
included as radiative corrections to the sum rules. Their inclusion is technically challenging, 
but does not pose a problem of principle. 

Our paper is organized as follows. In Sec. 2 we introduce the two different types 
of QCD sum rules. In Sec. 3 we then analyse and explain the discrepancy in B ^ p 
transitions, paying special attention to the rearrangement of quark and mixed condensate 
contributions within the light-cone expansion. We also estimate higher-twist contributions 
to Ai. In Sec. 4 we discuss shortly the quark mass dependence and the heavy quark 
limit. Section 5 contains our final predictions for form factors, spectra and decay rates of 
B —>■ pev, Sec. 6 a summary and the discussion of the results. More technical issues are 
delegated to appendices. 


^ While the usual sum rules are based on matching the QCD calculation at small distances to the 
phenomenological description in terms of hadrons at large distances, the light-cone sum rules match in 
transverse distance. Hence the relevant parameter in the expansion is twist, and not dimension. 
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2 Three-Point vs. Light-Cone Sum Rules 

In this section we discuss the two different types of QCD sum rules that have been used 
in the literature to calculate heavy-to-light meson decays. 


2.1 Three-Point Sum Rules 

The - comparatively speaking - “traditional” approach towards transition matrix elements 
is by calculating a three-point (“3pt” in the following) correlation function. Specihcally, 
for i? —p one consider^ 

ri*- = s"! {0\Tj;(v){V - y4)'‘(0)it (i)|0> 


= + pj'‘p^r+ - 


V 


+ 


( 2 , 1 ) 


Here (V—= -U7^(l—75)6 is the weak current mediating the b ^ u transition, js = qi'ybb 
is the interpolating held for the pseudoscalar B meson, and jp = the interpolating 

held for the p meson. In (|2.1| ) we have given explicitly only those Lorentz structures that 
are relevant for the semileptonic decay channel, the others are suppressed. 

The method of QCD sum rules consists - in principle - in performing on the one hand 
a QCD calculation of in the not so deep Euclidean region ~0(lGeV2), 

writing it on the other hand as (double) dispersion relation over the physical cut and 


equating both expressions. It was the idea of Shifman, Vainshtein and Zakharov to 
complement the purely perturbative calculation of T^^ by nonperturbative terms in form of 
vacuum matrix elements of gauge-invariant local operators, the condensates. The method 
proved surprisingly successful in describing various hadronic matrix elements in terms of a 
handful of input parameters, as is testihed by the immense number of publications in the 
held. 

The traditional approach by SVZ appeals to Wilson’s operator product expansion 
(OPE), which is the expansion of a T-product of currents at short distances in terms 
of local operators. In that way one obtains for the invariants T in (p.l|) : 


r(p|,pj. t) = 01 10) 


( 2 . 2 ) 


with the condensates ( 0 | | 0 ). In most applications one takes into account condensates 

up to dimension 6. T can also be expressed as a dispersion relation over the physical 
singularities: 

r= [dsBdsp—^ - f ’ ^ -l- subtractions. (2.3) 

J {sb-Pb){sp-p;) 

Usually one is interested only in the properties of the ground state, which has to be 
extracted from the sum over all states. For this, one writes 


^phys ^ground state i ^cont 

P = P + p 


^For 3pt functions we follow the notations of Q. 


(2.4) 
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and approximates the unknown contribution of the continuum by the perturbative spectral 
function above some “continuum-thresholds” Sq and Sq, such that 

^cont ^ (2.5) 

where the indicates that smearing over a sufficiently large interval is implied. The 
sum rule is then obtained by equating and ( p.3|) and subtracting from both sides the 
continuum contribution, i.e. in the above approximation the integral over the perturbative 
double spectral function above thresholds. In order to reduce the impact of this approxi¬ 
mation on the hnal results as well as the error induced by truncating the series (|2.2|) after 
the hrst few terms, one subjects the sum rule to a Borel transformation. For an arbitrary 
function of Euclidean momentum, f{P‘^) with the transformation is dehned by: 

(2.6) 

P'^/N=M‘^ fixed 


where is the so-called Borel parameter. For a typical term in the OPE, the transfor¬ 
mation yields 


Bp2{M^) 


1 

(p2 _ 


{n — 1)! 


■ir 


^—rrPjM'^ 


{M' 


2\n 


(2.7) 


As condensates with large dimension get divided by correspondingly high powers of — 
m^), their contributions to the sum rules get suppressed by factorials. In the dispersion 
integrals, gets exponentially suppressed relatively to the ground state contribution, 
which is just the desired effect. 

Dehning the couplings of the mesons to their interpolating helds by 


{0\dij5b\B^) = fB 


m 


B 


mb 


{0 I I P+, A ) = 


( 2 . 8 ) 


we hnd the following sum rules for the form factors determining the semileptonic B ^ p 
transition: 




mb 


fBfp{rnB + mp)m%mf 


exp 


ml 


m. 


M2 J ^ P 


mb{mB + rup) 

2 P = —FT—2-^ exp 

iBlpmlmp 




1/ 


^ mbjmB + rUp) 


2fBfpmlm^ 


exp 


ml ml 

M^^p 


ml ml 


MlM^ptp, 


MlM^pTv, 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


where the f on the right-hand sides denote the Borel-transformed invariants after contin¬ 
uum subtraction. The explicit formulas for F can be found in 


16 
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2.2 Light-Cone Sum Rules 


An alternative approach [^, 
the vacuum and the p meson state: 

Wp„A)|r(V'-yiv(o)il,(i)|o> 


2 ^ starts from the two-point function sandwiched between 




= + ^n2(pB,t)(e*(^Vs)(Ps +Pp)/A + ny(p|,f)e/“^e;:(^WPp/3 +(2-1?) 


with — ml < 0 and pi = m?. Hence we encounter only single dispersion relations, 


n = jdss ++£5?) + subtractions, 

J sb-Pb 


(2.13) 


and to isolate the ground state only an approximation for the continuum contribution in 
the B meson channel is needed: 


P 


cont 


~ pP«rt{i _ 0^^ 


0,(.B 


Sb)}- 


(2.14) 


Thus, this part becomes simpler and potentially more accurate than with 3pt sum rules, 
since less assumptions are made. 

The price to pay, however, is that the QCD calculation becomes more involved. In 


particular, the expansion of (|2.13|) in local operators becomes useless since an inhnite 
sequence of such operators contributes to the same order in l/(ps ~ bt^I). Indeed, each 
operator of the sequence 

g-Dpi ■ ■ ■-Dp^Tg, 

where is the covariant derivative and T an arbitrary Dirac matrix structure, sym¬ 
metrized over the indices pi,..., pn and with subtracted traces, enters with the same sup¬ 
pression factor I/{p's — Bxl) |^|. This is different from 3pt sum rules where contributions 


with higher n are supressed by powers of 1/Pp, which here is no longer an expansion param¬ 
eter. Still, contributions of operators containing traces over Lorentz indices, or transverse 
components of the gluon helds are suppressed by extra powers of l/(ps ~Brl). This means 
that the relevant parameter is the operator twist rather than dimension. The expansion 
goes in terms of nonlocal string-like operators on the light-cone (“LC” in the following), 
whose vacuum-to-meson transition amplitudes dehne the meson LC distribution ampli¬ 
tudes, which describe the momentum fraction distribution among the meson constituents. 
The leading-twist distributions correspond to the minimum number of Fock constituents 


and in the case of a charged p meson involve the following functions |18i : 


{p+{p,\)\u{0)a^^d{x)\0) = 


Pp) / dMe*“P>_L(M,/i), (2.15) 

Jo 


(p+(p, A)|M(0)7pd(a;)|0) = p^ 


(e++a;) 

{px) 


fpmp dne*“^’"^0||(M,/i) 
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{p-^{p,\)\u{0)jf,-f5d{x)\0) = -ef^^p^e 

For the sum rules for the B —> 


*ix)upp^af^mp [ due^'^P^g^^\u,p) 


1 

4 ' ■ '' Jo 

p form factors, we also need the function 

1 


(2.17) 


(2.18) 


In the above dehnitions the matrix elements are not gauge-invariant, but refer to the axial 
gauge XpA^{x) = 0. In a general gauge, gauge factors 

7-1 


Pexp 


^9 


da x^Ap{ax) 


(2.19) 


have to be inserted in between the quark helds. The integration variable u corresponds to 
the momentum fraction carried by the quark. The normalization is such that 

/ du f{u) = 1 
Jo 

for all four distributions / = (j)±, //(h); the scale-dependent coupling of the 

p meson to the tensor current, is dehned by (|2.15|) for x = 0. To leading-twist 2 accuracy 


the “^f-functions” are in fact not independent, but related to (f)\\ by Wandzura-Wilczek 
type relations: 


91 [u,p) = 


(a),twist 2/ \ o 

gy (M,/i) = 2 


dv 




+ 


dv 


h{v,p) 


u 


dv 


0||(u,/i) 


+ u dv 


011 (U,p) 


( 2 . 20 ) 


All these functions are discussed in detail in |^. The distribution amplitudes of higher- 
twist can be related to contributions of Fock-states with more constituents (say, an extra 
gluon) and generally generate contributions to the sum rules that are suppressed by powers 
of 1/{pb ~^b)- We will discuss higher-twist distribution amplitudes only shortly in Sec. 3. 

Note that the matrix element of a single nonlocal operator contains information about a 
whole series of matrix elements of local operators of increasing dimension (but hxed twist), 
which are encoded in moments of the distribution amplitudes. For example: 

{p-^ip, X)\u{0)apy{D ■ xyd{0)\0) = - el^^Pp)fp{ipxy f duu^(j)yu, p). ( 2 . 21 ) 

0 

A renormalization group analysis shows that for large n, moments of the above dehned 
distribution amplitudes behave as^ 

/ /i) ~ const/n^, (2.22) 

Jo 


Note that a purely perturbative analysis is sufficient to obtain the leading behaviour in n, whereas 
the coefficient of proportionality can only be obtained by using nonperturbative methods, see [03. 


































which corresponds to the following end-point behavior of the amplitude for n —> 1: 


const ■ (1 — m). 


(2.23) 


We would like to stress that it is in this place - using the information about large n 
behaviour of local operator contributions, related to the end-point behaviour of the LC 
distributions - that the LC sum rules go beyond the traditional SVZ approach. We will 
discuss this point in detail in the next section. 

The rest of the LC sum rule procedure follows the standard rules sketched in the last 
subsection: to suppress contributions of higher-twist p meson distributions and to enhance 
the sensitivity to the ground B meson state, one performs a Borel transformation in 
and the hnal expressions for the sum rules for i? —> p decay form factors are similar to the 
3pt sum rules in (|2.11|) apart from the different expressions on the right-hand side. One 
thus obtains (to leading-twist accuracy): 


A^r\t) = 


rrib 


fsimB + mp)mB 


2 exp 


. Ml 


B du 
/ — exp 
/o u 


u 


\uMl 


{t — ml — uml) 


0 [ c ( m , s ^)] \ fl-{p)(j)±{u,p) + + fp^bmpQ^l^u, p) [>, 


(2.24) 


mb{mB + mp) f 

A if) = - 7 - 2 - \ 


m%-ml 

Ml 


B du 

/ — exp 
10 u 


u 


\uMl 


{t — ml — uml) 


fp{f(p±{u, p)Q[c{u, So )] + fpmbmp^w {u, p) 


uMl 


0 [ c ( m , s ^)] +6[c{u,sf] 


>, (2.25) 




mb{mB + mp) j' it^b ~'fb\ f du j 


u 


2fBm 


exp 


B 


Mi 


— exp 
0 u 


[uMl 


(t — ml — umf) 


fpif(t>fu,fQ[ciu,sf] + l^fpmbmpg^l\u,p) 


0 [ c ( m , s ^)] + 5 [ c ( m , s ")] 


B\ 


uMl 


>.(2.26) 


with c ( m , sf = uSq — ml + tu — uuml- 

LC sum rules for Ai and V were already obtained in |^; they slightly differ from the 
ones given above by the “surface” terms ( 5 [ c ( m , sU which are related to subtleties in the 
continuum subtraction as discussed in App. B. The LC sum rule for A 2 is new. 


2.3 The Conflict 

The two approaches described above are rather different and their comparison should shed 
light on the actual accuracy of the sum rule method. The numerical comparison requires the 
use of a “coherent” set of parameters, so that differences are not introduced (or masked) 
by using different inputs. We shall specify our set in detail below; for the purpose of 
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0 5 10 15 20 0 0.2 0.4 0.6 0.8 1 1.2 

r[GeV^] r[GeV^] 



Figure 2: Semileptonic form factors of the decays D ^ p and B ^ p from LC sum rules, 
Eqs. (|2.24|) - (|2.26|) , and 3pt sum rules, Eqs. (|2.9|) - (|2.11|) , evaluated with the same input 
parameters. The dotted curves illustrate the effect of introducing a different interpolating 
current for the p meson in the 3pt correlation function, see Sec. 3.3. 
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illustration the particular values are unimportant. The results for all B and D meson 
decay form factors from both 3pt and LC sum rules are shown in Fig. 2. We see that the 
results are in reasonable agreement at large t while there is a disturbing discrepancy up to 
a factor two at large recoil [^, |^. The t dependence also turns out to be very different 

i, Un¬ 


provided no particular advantage or flaw of one method can be found, this spread 
of values would necessarily have to be considered as indicating poor theoretical accuracy 
of the predictions in this region. The further discussion will clarify the reason for this 
discrepancy and give strong evidence in favour of the LC sum rule calculation. Reasons 
for the better agreement at small recoil (large t) will also become clear. 


3 Anatomy of the Discrepancy 

An inspection shows that the disagreement between LC and 3pt sum rules is mainly due 
to the contribution of the quark condensate, which dominates the 3pt sum rules at small 
t, cf. [^, 1^. To clarify the reason, we give in this section a detailed calculation of this 
contribution, and also of the contribution of the mixed condensate to the 3pt sum rule for 
the axial form factor Ai. The result is well known ||T^ , ^ and the new point we wish to make 
here is to rederive it using the sequence of steps adopted by the LC sum rule approach. 
This will reveal how the p meson distribution amplitudes are implicitly described in the 
3pt approach and also give examples of higher-twist contributions. 


3.1 Three-Point Sum Rule from the Light-Cone Point of View 


We start from the correlation function (^) and as hrst step substitute the heavy quark 
propagator by its leading-order perturbative expression: 


d^x d^y J 


d*k 

{ 271^1 


Akx 


ml — k‘^ 


X {0|rR(9)u(0)7'‘(l -75)(mi,+ ^)s75(i(a:)} |0). 


(3,1) 


The product of q-matrices in (p.l|) contains several terms, corresponding to different in¬ 
variant structures in (^) and to contributions of dimension odd (even) operators to the 
OPE. We choose to consider the axial form factor Ai, and contributions of operators of 
odd dimension only. To this end we need to calculate the correlation functions 


, I dSe^py (0|T{j;(2/)u(0)a^“d(a;)}|0), (3.2) 

= ^| d%e^Py {0\T{j;{y)u{0)d{x)}\0), (3.3) 

using the OPE (we recall that is assumed to be Euclidian and sufficiently large). 
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Throughout the calculation we imply using Fock-Schwinger gauge. In a general (co¬ 
variant) gauge the heavy quark propagator in external gluon helds contains the link factor 
( p.l9| ), which has to be inserted in the nonlocal operators in (|3.2|) and ( |3.3| ) to make them 
gauge-invariant, see Sec. 2.2. In the Fock-Schwinger gauge, further terms in the expansion 
of the b quark propagator in background gluon helds only yield corrections ~ l/mjy to the 
sum rules and for simplicity will not be considered here. They can easily be added. 

The OPE of the correlation functions ( p.2|) and (|3.3|) is straightforward and yields:^ 




r 1 


1 ml 


1 1 


^ — 




3 p 


16 






1 

px 


+ t{qq)—^{x^gP^ 


1 7T7 ^ 
6 p^ 


1 T)^ 1 777 ^ 

+ i{qq) — ^{p'^x^ - p^x'^){ipx){l - , 

px p^ 6 p^ 


(3.4) 


rj^U _ 


i{qq) {p''{px) — x''p^^ (1 -|- 


6 p 




p. 


2 2 \ 
m^x 1 . 


(3.5) 


Here = {qgcGq)/{qq). Note that = 0, while contains a contact term. 

Substituting o and (|3.5|) into (0. taking the remaining integrals and performing 
Borel transformations in and p^, respectively, we reproduce the contributions of quark 
and mixed condensate to the 3pt sum rule for Ai in Ref. [Q, except for the neglected 
contribution of the diagram with the gluon emitted from the b quark line: 


Hr^(t) = 


mb 


m 


exp 


B 


mt 


Ml 


fnlpirriB + mp)m%mp 

mliml — t) ml — t 


+ 


{qgaGq) 


7 


Ml 


3ml 


-{qq)^iml-t) 


8Ml 


QM^ 


24M| 


5t ^ (m\ 


tf 


6MlM^ 


(3.6) 


We emphasize that the derivation sketched above is entirely within the traditional QCD 
sum rule approach, although the sequence of steps may seem unusual. A related discussion 
for the 7r77* transition form factor can be found in [TO . 


We now rewrite this answer in terms of contributions of p meson distribution ampli¬ 
tudes. To this end, we separate the p meson contribution to T^^"(p), 


= m\p^{pA))- 


mt 


pz 


:{p+(p,A)|ti(0)CT'‘“d(i)|0> + .... 


(3.7) 


and, similarly, the one to T^. The first matrix element is proportional to the decay constant 
fp, while the second one, by definition, gives p meson distribntion amplitudes in the fraction 


The perturbative contribution to (3.2) and (3.3) is of order mu,d and will be neglected here. 
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of momentum carried by the quark. An inspection of ( |3.4| ) and ( pi.51) suggests to introduce 
the following distributions: 

(p+(p, A)|h(0)cT^„d(a;)|0) = - e*J^'>Pp) j due^'^P^[(j)^{u) + 

0 

+ ■ x){xpPa-XaPp) [ due^'^P^'ijj^^\u), (3.8) 


(^■'■(p, A)|u(0)(i(a;)|0) = —■ x) [ du 


(3.9) 


After a Borel transformation of (|3.4|) and ( p.7|) in p^ we get the explicit expressions 
4>x{u) = 


■im) ^rallMl if, 1 "^0 


■ sifj ^ - “4 

(3.10) 


y‘’(“) = —A™o|i5(“) +i5(1 

^pjpjp J-t* 


= —7^e”’?''"’lm5[i5(u) - i5(l - ti)], 

’"p/c/o 6 


(3.11) 

(3.12) 

(3.13) 


'ip^‘^\u) = 0 , 

where (1 — 2) GeV^ is the Borel parameter. Note that the expansion goes in deriva¬ 

tives of the 5-function. A similar expansion for the twist 2 distribution amplitude was 
obtained in . 


Similarly, from the expansion (|3.5|) we deduce 




{qq) ^2/,^2l 


^pfpfp 


- ml[6{u) + 5(1 — u)]. 
6 


(3.14) 


Substitjuijin^ Oiiici iri Sinci (EH). taking the integrals and performing a 
Borel transformation in p^, we get a typical LC sum rule: 

rribf^ „„„ f 

\ o 6Xp ^ ) j 

JO u 


A^r^t) = 




fB{mB + mp)m% 

0 u(m 


u 


Ml 


exp <i -772 UTTip) 


[uMl 


uMl 


1 + 


uMl^ 


'ip^^\u) H- 772 '^^‘^\'^) 


uM 


B 


+ 


U 


1 + 


2ml 

uMl 




1 + 


2ml 

uMl 


u 
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(3.15) 

























where we have changed variables m —>■ 1 — m to be consistent with ( ^.24 ). To save space 
we have not shown the continnnm snbtraction. Note that the leading-twist contribntion 
of the distribntion amplitnde (j)± coincides with the corresponding contribntion in ( p.24| ); 
the extra terms "0^*^ are higher-twist corrections, not taken into acconnt in (p.24|) .f] 


On the other hand, fnrther snbstitnting in (|3.15|) the above expressions for the distri¬ 


bntion amplitndes and snppressing terms ~ nip nnless they get divided by we come 
back to (|3.6|) . The qnark condensate contribntion in (^T|) appears as a contribntion of 
the leading-twist distribntion (j)±, while the mixed condensate terms contain contribntions 
from both leading- and higher-twist. In particnlar, for the expression in sqnare brackets in 
(p.6|) we hnd the following decomposition: 


{ml — ty ml — t 


QMlM^p 




+ 


M^p 


6 + 


ml 


— t 


^Ml 


+ 


-I 


ml{ml — t) ^ ml — t 


8Mt 


B 


8M^ 


B 


4M| 




ml 

3M| 




'0(4) 


(3.16) 


where [.. indicates that this term originates from the distribntion 0^, etc. As it stands, 
this expression does not yet agree with ( p.6|) , the reason being that the Borel transformation 
in the p meson channel was applied in a slightly different way. It is possible to show that 
in order to reprodnce the 3pt snm rnie one has to snbstitnte —*■ —Mp, after which the 
expressions indeed coincide literally^. 


3.2 Short Distance Expansion and Light-Cone Distribution Am¬ 
plitudes 

The new inpnt made by the LC snm rnles is to argne that the (5-fnnction type shape of 
LC distribntions, concentrated at n = 0 and n = 1, is qualitatively wrong. In particnlar. 


instead of the expression in (|3.10|) , it is snggested to nse the distribntion amplitnde 


0 ±('W, p) = 6n(l — u) 1 -|- a 2 {p)C 2 ^‘^{2u — 1) 


(3.17) 


Here Cl^‘^{x) 


= (15a;^ ~ 3)/2 is the second order Gegenbaner polynomial; the coefficient 


was estimated to be 0.2 ± 0.1 [0. Eq. (|3.17|) is clearly very different from ( p.lOj) . Where 


does it come from and what is wrong with (p.l0|) ? 

The distribntions (|3.10|) - (p.l4 ) are jnst the QCD snm rnles for the correlation fnnctions 
(p.2|) and (p.3|). Their dehciency becomes apparent when they are rewritten in terms of 


® The contribution ^ g± in ( p.24| ) would correspond to terms in S) with an odd number of y-matrices, 
which we have not considered here. 

^ There is a subtlety in treating the terms proportional to Pajp^ in the hrst line in (3.4): pa gets 


contracted with ps and yields a factor p^^ -\- p'^ — t. Using the dispersion relation first in the p meson 
channel like in (0) then implies that is substituted by rUp, while in the standard procedure it gives 
—Mp. Ambiguities of this type are intrinsic for the sum rule method. 
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Figure 3: Ai from the LC sum rule ( p.24| ) (solid line LC), the 3pt sum rule ( p.9|) (solid line 
3pt), and a “hybrid LC” sum rule, in which the leading-twist distribution amplitude (f)± is 
replaced by the naive expansion in ^-functions, (|3.10|) . 


moments. For the leading-twist distribution we hnd (cf. 


[ du (2u — l)'^4>±(u) = 
Jo 


-(gg) 

^pfpfp 


jmllMl 


(i + (-in 1 


^5 

3M2 


(2n + 1) 


(3.18) 


Note that the contribution of the mixed condensate is enhanced by a factor n. This 
enhancement is of generic nature: contributions of higher dimension to the OPE will be 
acompanied by increasing powers of n so that the sum rule blows up for large moments 
and cannot be used. This signals the break-down of local OPE for higher moments of 
distribution amplitudes. Extensive studies ||2^ have demonstrated that QCD sum rules of 
type ( |3.10|) - (|3.14|) can be applied to estimate the two hrst moments only, n = 0 and n = 2, 
i.e. the normalization and width of the distribution amplitudes, but fail to describe higher 
moments, i.e. the shape of the distribution close to the end-points. Information on the 
shape can, however, be obtained from another source, namely the behaviour of distribution 
amplitudes under the renormalization group [O. The major result is that (j)± approaches 


6 n(l — m) at large virtualities and that the corrections can be systematically expanded in 
Gegenbauer polynomials C^^‘^{2u — 1). Combining this expansion with estimates of the 
hrst two moments by QCD sum rules one obtains the expression (|3.17|) . 

In fact, the particular sum rule in (|3.18|) is not accurate enough even for n = 0, 2, and 
in practice one uses different sum rules, see Ref. 


18 for a detailed discussion. 


To illustrate that the shape of the leading-twist distribution is indeed of crucial impor¬ 
tance, we have plotted in Fig. 3 the form factor 74f^^(t), calculated in several different 
ways. The solid curve, labelled LC, shows the LC sum rule (p.24|) with realistic distribu¬ 
tion amplitudes. The dotted line is obtained using the same sum rule (^.24), but with the 


distribution amplitude (j)± replaced by the expression (|3.10|) ; it is very close to the solid 
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line showing the 3pt sum rule result. The “dominance” of the quark condensate in the 


3pt sum rule thus happens to be an artifact of the short-distance expansion extrapolated 
beyond the region of its validity. 

The ideal agreement of the dotted curve in Fig. 3 with the 3pt sum rule result at f = 0 is 
in fact coincidental and is due to a mutual cancellation of two effects. First, in addition to 
contributions of operators of odd dimension, the 3pt sum rule contains a perturbative term, 
a contribution of four-quark operators of dimension 6, and of the gluon condensate. These 
contributions correspond to the terms with an odd number of 7 -matrices in (|3.1|) , which we 
have not considered, and have their counterpart in the LC sum rule in the contribution of 
the distribution g± (up to higher-twist terms). The difference between the two approaches 
is small in this case, the reason being that repeating the above procedure one would deal 
with the correlation function of with a nonlocal vector current. In contrast to ( |3.2| ), 


( p.3|) it has a large perturbative contribution and the OPE goes in condensates of even 
dimension. Extracting the distribution amplitude as outlined above would yield a smooth 
distribution ~ u{l — u), slightly corrected by ^-function type contributions of the gluon and 
four-quark condensates. These latter contributions are small, so that g± as implicitly used 
in the 3pt sum rules is not very different from its “true” behaviour. Hence, the numerical 
results are close. 

Second, the present version of the LC sum rule neglects contributions of higher-twist. 
To estimate their effect one can apply the methods of Ref. to determine the shape of 
the distributions at large scales, i.e. their asymptotic form, and use the sum rules 

( P.13| ) to estimate the normalization. We get 

= Ki ■ 30n^(l — m)^, 

= K 2 ■ 30n(l — m)(1 — 2m), 


= K4 ■ 6m(1 — u) 


with 


Ki = ml/16, K 2 = rnl/G, H 4 = —ml/ 6 . 


(3.19) 


(3.20) 


In Fig. 4 we plot from Eq. (|3.15|) using these distributions and including continuum 


subtraction. For comparison we also show the leading-twist LC sum rule ( p.24|) . The 
correction turns out to be negative and lowers the leading-twist result by about 15% for t < 
15 GeV^. These results are, however, only indicative on the size of higher-twist corrections, 
the detailed study of which goes beyond the tasks of this paper. 

If the “naive” description of distribution amplitudes by the usual sum rule method is 
that dehcient, the question arises if this approach still works for form factors of D mesons, 
as used e.g. in |]I3. The formal answer is clear from the structure of LC sum rules: the 


distribution amplitudes are integrated with a smooth weight function over a constrained 
region of the momentum fraction u. If the mass of the heavy meson is not very large 
compared to the typical hadronic scale 1 GeV, then the integration region is large and only 
gross characteristics of the distribution amplitudes matter, i.e. their normalization and 
width. These are given correctly by the sum rules, and the approach works well. If, on the 
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Figure 4: Ai from the LC sum rules (|2.24|) and (|3.15|) , with the distribution amplitudes 


other hand, the mass of the heavy meson is much larger than 1 GeV, as it is the case with B 
mesons, and if the momentum transfer to the leptons is small, then the integration region 
shrinks to the narrow interval 1 — u ~ 0 {l/mb), the precise behaviour of the distribution 
amplitude at u —> 1 becomes important, and the standard approach fails. 

The physical parameter that matters is, however, not the heavy meson mass, but the p 
meson energy Ep in the decaying B{D) rest frame: Ep = — t)/{2mB,D)- Zero 

recoil corresponds to Ep = nip-, in the physical regime t > 0, Ep runs up to 2.7 GeV and 
1.1 GeV in B and D decays, respectively. In Fig. 5 we show the form factors Ai{Ep) for 
both B and D mesons. The behaviour is very similar, and in both channels 3pt and LG 
sum rules agree very well for Ep 1.4 GeV. For D mesons, this is outside the physical 
region for the decay. 

3.3 Possible Remedy: the Tensor-Current? 

To conclude this section, we would like to demonstrate that the “dominance” of the quark 
condensate is no intrinsic feature of 3pt sum rules. To this end we recall that one has 
some freedom in the choice of the interpolating held for the considered particles: although 
for the p meson the vector current is the most convenient one, it is by no means the only 
one. In particular one can choose the tensor current instead and calculate the 

B ^ p form factors from the correlation function 

=fj {0\Tfp^{y){V - 7l)^(0)j|j(x)|0). (3.21) 

In App. B we give the corresponding OPE including terms up to dimension six. Due to 
the particular y-matrix structure, the contribution of the quark condensate to A2 and V 
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Figure 5: Ai as function of the p meson energy Ep from LC and 3pt sum rules for both 
D ^ p and B ^ p transitions. For D ^ p transitions, the physical region t > 0 
corresponds to p meson energy up to 1.1 GeV. 

vanishes and is small for Ai. We have displayed the corresponding form factors already in 
Fig. 2. They differ distinctively from the results of the original 3pt sum rules and are much 
closer numerically to the LC sum rules. Nevertheless it would be inappropriate to conclude 
that the above correlation function is “better” than (|2.1|) : it suffers from exactly the same 
problem as the original correlation function to describe correctly the shape of the p meson 
distribution amplitudes near the end-points in u. It is only that this failure is less “visible” 
for the given values of the b quark mass and the considered range in t. The problem is 
now shifted to the contribution of the mixed condensate, which starts to domininate ^i(O) 
at large and eventually overgrows all other terms. Numerically, however, the effect is 
much less signihcant at m;, ~ 5 GeV. This improvement comes at the price that the tensor 
current couples also to positive parity states which contaminate the contribution of 
the p meson, so that the accuracy of these sum rules is not very high. Another possibility 
to achieve a similarly “favourable” rearrangement of power corrections would be to use the 
axial-vector instead of the pseudo-scalar current for the B meson. 


4 The Heavy Quark Limit 

The behavior of form factors in the limit mi, —cx) is of considerable theoretical and 
practical interest. Taking the heavy quark limit in the sum rules is straightforward, by 
rescaling the sum rule parameters in the following way (see e.g. |3^ ): 

^ 2mfer, ml + 2mbUJo, (4.1) 

where r and Uq are of order 1 GeV. 
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One should distinguish between different regions of momentum transfer. First, consider 
ml — t ~ i.e. small recoil, energy of the outgoing p meson of order 1 GeV. Then 

both the 3pt and the LC sum rules fulfil the scaling laws predicted by heavy quark effective 
theory 


Ai{t ~ 

tmax ) ~ A2{t 

tynax ) ~ v{t^ tynax ) 


(4.2) 


In this regime, the integration over the quark momentum fraction in LC sum rules comprises 
the region 1 — u ~ 0 ( 1 ), so that only width and normalization of the distributions are 
important. Hence 3pt and LC sum rules are expected and indeed give comparable results, 
see Figs. 2,3. 

More interesting, however, is the behaviour near maximum recoil, t ^ 0. Here we find 
that in the 3pt sum rules approach the limit mb oo cannot be taken since higher order 
terms in the OPE are accompanied by increasing powers of mb. From the “light-cone point 
of view” this inconsistency arises because at large recoil the soft contributions to the form 
factors pick up a tiny region of momentum fraction 1 — m ~ 0 (l/m;,) and thus the details 
of the shape of the p meson distribution amplitudes, wrongly described by 3pt sum rules, 
enter decisively. 


On the contrary, LC sum rules at t = 0 have a well-defined heavy quark limit and 
scale as Explicitly, making the change of variables a; = (1 — u)mb/2, one finds 

(with / = /syfub and A = ms — mb). 


Mi(0) = 

fA2{0) = 

fV{0) = 


m, 


3/2 

b 

2 


m, 


3/2 

b 

2 


^A/r 


^A/r 


fU>0 


ru>o 


ru>o 


duje - fpmpg^l\l) 

due--/- [/W(l) + /,m,$j|(l)^ 


m, 


3/2 


/ due-/- 


+ ^fp'mp9T\^) 


From the relations (|2.18|) and ( p.20| ) it follows that to our accuracy 

Hi(0) = H2(0) = H(0) 


(4.3) 


(4.4) 


in the heavy quark limit. This agrees with the findings of Ref. |^. It is instructive to 
check that the above scaling relations are not spoiled by higher-twist corrections. The twist 
4 part of (|3.15|) becomes in the heavy quark limit; 


/H™(0) = 


C ‘^0 


m 


5/2 


/ du e-l- 


.i)/!v,"4)(l) - , (4,5) 


r 


r 


where 'ip'{u) = {d/du)'ip{u) and we used that all //'-functions vanish at m = 1 . It is seen 
that higher-twist corrections are in fact down by an extra power of mb, cf. the discussion 
of the pion form factor in the third of Refs. . 
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We recall that the heavy quark mass dependence of form factors at zero recoil is of 
vivid interest for lattice calculations. Due to restrictions on computer power and perfor¬ 
mance, reliably simulable quark masses are of order ~ 2 GeV and the results have to be 
extrapolated to the physical beauty quark mass. In this respect, we would like to add 
a word of caution about using the asymptotic scaling law since this limit is only 


approached very slowly [O. To get a ball-park estimate of the next-to-leading order cor¬ 


rections we calculated the form factors using LC sum rules varying the b quark mass in 
the limits (1-10) GeV and using the scaling ([4.2|) of the sum rule parameters. We then £t 
by a quadratic polynomial in the inverse meson mass ms = mt, -|- A, A = 500 MeV ^ . 
The results are (we show the leading l/niB corrections only): 

mfMO) = 5.6 (1 - + ,,, 

V niB 


m 


3/2 

B 


A2{0) = 5.6GeV^/2 (l - 


2.1 GeV 


niB 


+ ... 


m|/V(0) = 5.9GeV^/2 (l 


1.5 GeV 
tub 


+ ... 


(4.6) 


The constants in front of the brackets are almost equal, as expected from ( |4.4|) . Note the 
large terms in l/m^. 

Finally, one can consider the region of very small recoil ml — f ~ 0(1 GeV). This 
region is generally difficult for QGD sum rule treatment since one gets more sensitive to 
contributions of large distances in the “t-channel”. An inspection of ( 3.16|) shows that in 
this limit the leading-twist contributions of dimension 5 are smaller than those of higher- 
twist, which may be considered as an indication that 3pt sum rules become more reliable 
than LG sum rules at very large t. 


5 Numerical Analysis 


We now turn to the numerical evaluation of the LG sum rules (|2.24|) -( |2.26|) . Let us first 
dehne the relevant observables. 


5.1 Kinematics 


With the standard decomposition for the B ^ p transition matrix element ( p..l|) the 
spectrum with respect to the electron energy E reads: 

dT{B^ p+e-jy) _ 
dE “ 


2 J dtt{{l — COS+ {1 + COS + 2{1 — cos^ 6 )Hq}, (5.1) 
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with the helicity amplitudes 


H± = {iriB + mp)Ai{t) =F 


Ai /2 


niB + rn 




Hn = 


2 mpy/t 


(m| - ml- t){mB + mp)Ai{t) 


A 


mB + mp 




where the indices denote the polarization of the p. A is dehned as 


A = (rn\ + 

^max, the maximum value of t at hxed electron energy, is given by 

.2 \ 


tynax [ mB 


m„ 


mB — 2E 


(5.2) 

(5.3) 


(5.4) 


(5.5) 


6 is the angle between the p and the charged lepton in the (e z/) CM system and given by 

1 


cosO = 


Ai /2 


{m^ — mi + t — AmBE). 


The spectrum with respect to t reads: 


dT{B^ 


p'^e u) 


Gj,\V, 


ub 


dt 


192TT^m% 


t (hI + hI + Hpj . 


(5.6) 


(5.7) 


We also introduce the notations Tt and T^ for the partial decay rates where the hnal state 
p is transversely or longitudinally polarized. 

From the specihc structure of the helicity amplitudes it follows that at small t the 
produced p mesons are predominantly longitudinally polarized; for f = 0 only longitudinally 
polarized p are produced. At large t, on the other hand, the contribution of A 2 and V to 
the decay rate is suppressed, since A has a zero at tmax- 


5.2 Input Parameters 

The decay constant fp is measured experimentally [P7|]: 


fp = (205 ± 10) MeV, 


(5.8) 


while existing information on comes from QCD sum rules. In the following we use |T^ 


f^{l GeV) = (160 ± 10) MeV. 


(5.9) 


The p meson leading twist distribution amplitudes and have been recently reexam¬ 
ined in 1^. We use 


0||,±(M,At) = 6 m(1 


u) 


1 a!^^{p)Cl^^{2u 


1 ) 


(5.10) 
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as already mentioned 


with a^(lGeV) = 0.18±0.10 |2|, 0 and a^(lGeV) = 0.2±0.1 ^ 
in Sec. 3.2. 

The value of the b quark (pole) mass mb is somewhat controversial, with estimates 
varying from 4.6 to 5.1 GeV. This large range, however, probably overestimates the actual 
uncertainty and rather reflects that the pole mass has to be nonperturbatively dehned and 
that suitable dehnitions (and values) depend on the application. In this paper we use 


mb = (4.8 ± 0.1) GeV, 


(5.11) 


which, we believe, is a fair estimate. 

The decay constant fs was calculated in QGD sum rules and on the lattice, with a 
world average of about 180 MeV (see, e.g. [^). It was found, however, that within the 


QGD sum rule approach receives large radiative corrections, which increase its value by 
30 to 60 MeV . Since similar radiative corrections have not been calculated for the sum 


rules for form factors, we think that it more consistent to use the lower value of fs as it 
is obtained from the sum rules without radiative corrections, see also 1^. In practice, we 


simply substitute by the corresponding sum rule with the same values of all parameters; 
this has an additional advantage of reducing considerably the b quark mass dependence. 
In fact, there are arguments suggesting that radiative corrections tend to cancel between 
Jb and the form factors. This cancellation was indeed observed for B D* transitions 
|3^ and for the B meson matrix element of the kinetic energy operator [Q. An explicit 
calculation of the radiative corrections to LG sum rules would, however, be very welcome. 

For the values of the condensates we use 


(gg)(l GeV) = -(245 ± 10) MeV^ 


TT 


= ( 0.012 ± 0.006) GeV^ 


{qgaGq){lGeV) = 0.65 GeV^ ■ (gg)(lGeV), 
{asqqf = 0.56 ■ (-0.245) GeV®. 


(5.12) 


They enter the 3pt sum rules explicitly, and the LG sum rules implicitly, via estimates of 
the parameters of the distribution amplitudes and of fs- 

We assume values of the continuum thresholds for p and B mesons Sq = 1.5 GeV and 
Sq = 35, 34, 33 GeV^ for mb = 4.7,4.8, 4.9 GeV, respectively. The working region of Borel 
parameters in 3pt sum rules is taken to be (IW) GeV^ for p mesons and (5- 

10) GeV^ for B mesons, with a hxed ratio M^/M^ = 5. Since for hxed momentum 
fraction u the expansion in LG sum rules goes in powers of 1 /(mM^), we make the formal 
replacement —>■ M^/ (u), where (u) ~ 0.6 — 0.8 is the average momentum fraction 

calculated by inserting an additional factor u under the integral (separately for each form 
factor and each value of t), and then taking the interval (4-8) GeV^, the same as in 

the 2pt sum rule for fs- The scale of condensates and distribution amplitudes in the sum 
rules for the form factors is = m\ — ml. 


22 
































5.3 Results and Error-Estimates 

Our final results for form factors and spectra are collected in Figs. 6 - 8 . First, we display 
in Fig. 6 the form factors as functions of the Borel parameter at f = 0. The solid lines are 
obtained with mb = 4.8 GeV {sq = 34GeV^), the dashed lines with mb = 4.7 GeV {sq = 
35 GeV^) and = 4.9 GeV {sq = 33 GeV^), respectively. The curves are remarkably flat 
which indicates a good accuracy of the sum rules. The variation of a|’^ within the specified 
±50% range has an effect of about the same size as the dependence on mb- The dependence 
on the continuum threshold Sq is small, provided the same value is used consistently in 
the sum rule for fs- In addition to uncertainties in the sum rule parameters, the accuracy 
of our results is essentially limited by the neglected higher twist corrections and radiative 
corrections. We have estimated the higher-twist effects for Ai in Sec. 3.2 and found them to 
be approximately —15%. This estimate is, however, preliminary and we have not included 
the higher-twist correction in our final results in this section. As for radiative corrections, 
we expect them to cancel to some extent when /b is expressed as 2 pt sum rule to 0 (a°) 
accuracy. Both sources of uncertainty can be systematically reduced by calculating the 
corresponding corrections, which is possible, but beyond the scope of this paper. Taking 
everything together, we think that adding an additional ±15% uncertainty to the above 
results yields a fair estimate of the true theoretical error. 

We thus obtain the following values for the form factors at maximum recoil: 

ylf^^(O) = 0.27 ± 0.01 ± 0.02 ± 0.02 ± 0.04, 

yl^^^(O) = 0.28 ± 0.01 ± 0.02 ± 0.02 ± 0.04, 

= 0.35 ± 0.01 ± 0.03 ± 0.03 ± 0.05, (5.13) 


where the first error comes from the variation in the Borel parameter, the second from the 
uncertainty ±0.1 GeV in mb, the third from the uncertainty ±0.1 in and the forth from 
the estimated uncertainty due to not included higher twist and radiative corrections. Note 
that the first three errors are correlated between the form factors. The results for Ai(0) 
and V(0) are comparable with those obtained in |]^. In Table 1 we compare our results 


to quark models, adding the errors in quadrature. We have not included the 3pt sum rule 
results 1 . 0 . since they suffer from the deficiencies discussed in Sec. 3. A comparison with 
lattice results is difficult, as most of them are obtained at large t > 14GeV^ and then 
extrapolated down to f = 0 using different assumptions on the functional dependence on t 
and the b quark mass. Only for Ai the assumed monopole dependence Ai ~ l/(m|* - t) 
is compatible with the scaling law Ai(0) 
collaborations have obtained 


Using that dependence, different lattice 


0.22 ±0.05 ELG [7], 
Ai(0) = { 0.24 ±0.12 APE g, 

UKQGD |n 


0 

^•^<- 0.04 


(6.14) 


These numbers are quite close to our result. 
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Figure 6 : Semileptonic B decay form factors at t = as function of the Borel parameter for 
central values of the LC sum rule parameters (solid lines). 
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Figure 7: Semileptonic B decay form factors as function of t for central values of the LC 
sum rule parameters (solid lines). The dashed lines give error estimates. 
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yB^P 

0.35 ±0.07 


Reference 
This work 


BKR 


FGM 41 


Jans 42 


Melikhov 
WSB 


11 


— >7r 


0.30 

0.20 ±0.02 
0.27 
0.29 
0.33 


AB^p 

0.27 ±0.05 

0.26 ±0.03 
0.26 

0.17-0.18 

0.28 


0.28 ±0.05 

0.31 ±0.03 

0.24 

0.155 

0.28 


0.29 ±0.03 
0.35 
0.215 
0.33 


Table 1: The form factors of the 6 —> m transitions at t = 0 in LC sum rules and quark 
models. 


Reference r(R0 ^ vr+e-p) r(R0 ^ p+e-p) r(p)/r(7r) Tl/Tt 


This work 

- 

13.5 ±4.0 

1.5± 0.5 

0.52 ±0.08 

BKR i 


8.7 

- 

- 

- 

FGM 

41 


3.0 ±0.6 

5.4 ± 1.2 

- 

0.5 ±0.3 

ISGW2 H 

9.6 

14.2 

1.48 

0.3 

Jans p2| 


10.0 

19.1 

1.91 

0.82 

Melikhov 

43 

7.2 

9.64 

1.34 

1.13 

WSB 

0 


7.4 

26 

3.5 

1.34 


Table 2: Decay rates of the b ^ u transitions in units |14bpps T^ denotes the portion 
of the rate with a longitudinal polarized p and Tt with a transversely polarized p. 


Next, we display in Fig. 7 the behaviour of the form factors in t (solid lines) together 
with error estimates (dashed lines) obtained by taking extreme values of the parameters: 
the upper dashed lines refer to = 4.7GeV, = 4GeV^, the lower dashed lines to 

rrih = 4.9 GeV, = 7 GeV^. We also show lattice results from the UKQGD collaboration 
(diamonds), which are in very good agreement with our results. The plots indicate clearly 
that the accuracy of our results at large t is worse than at small t, in particular for A2 and 
V. However, the contribution of A2 and V to the experimentally measurable observables, 
the spectrum in t e.g., is kinematically suppressed at large t, so that large uncertainties in 
that region are not relevant phenomenologically (see also the discussion below). Figure 7 
also shows that Ai is a slowly varying function of t, whereas A2 and V increase more 
steeply; none of the form factors can be fitted by a monopole in t as suggested by the 
pole-dominance hypothesis. In Ref. it was found that the ratio of form factors takes a 
simple form in the heavy quark limit supplemented by some model-assumptions. We find 
that in the full range of physical t our ratio V{t)/Ai{t) agrees with the prediction of 
within 4%, whereas A 2 (t)/Ai(t) is by 10% to 20% smaller than predicted. 

Finally, in Fig. 8, we show the spectra dV/dt and dV/dEe- Fig 8(a) shows the effect 
mentioned before: although the uncertainty in the form factors increases with t, the con¬ 
tribution of A2 and V is suppressed and the resulting uncertainty is dominated by the 
(smaller) error on Ai. The uncertainty is maximal at t ~ 15 GeV^ and amounts to 1^5%, 
which yields a accuracy of |I4b| if determined from that point. Taking into account 
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Figure 8: i? —>• pev decay spectra, (a): spectrum in t, (b): spectrum in the electron energy 
Ee- Solid and dashed lines as in Fig. 7. 


the additional uncertainties of unknown higher twist and radiative corrections, we estimate 
that with present knowledge \Vub\ may be determined from dT/dt with a theoretical accu¬ 
racy of 20%. It is conceivable that further calculations may push down this uncertainty to 
±15% on the spectrum, i.e. about 8% on \Vub\, especially if mb was fixed to better accuracy. 
Fig. 8(b) also shows that a determination of \Vub\ from the electron energy spectrum may 
be more difficult, since it is strongly peaked and the position of the maximum thus may 
be invisible with presently available experimental resolution. 

In Ref. PI the CLEO collaboration has presented first results on the branching ratios of 
B —>■ Tieu and B —>■ pev. Since the given values are to a certain extent model-dependent, 
we refrain from extracting any number for \Vub\ from them. This task, we believe, is more 
appropriate for our experimental colleagues. 

Integrating up the spectra, we hnd 

F(R° ^ p+e-p) = |14feP(13.5 ± 1.0 ± 1.3 ± 0.6 ± 3.6) ps"^ (5.15) 

with the same sequence of errors as for the form factors. In Table 2 we also give ratios of 
partial decay rates which are independent of \ Vub\ and may serve as tests of our predictions. 
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To get the ratio r(p)/r(7r) we have used the result of [Q obtained by a similar method. 


6 Summary and Conclusions 


We have given a detailed analysis of existing controversies in QCD sum rule calculations 
of semileptonic B pev form factors, which, as we believe, settles this problem. Both the 
decease of 3pt sum rules, which we have exposed, and the remedies which we have suggested, 
apply to all heavy-to-light transitons and are equally relevant e.g. for rare radiative decays, 
where a similar discrepancy between LC and 3pt sum rules was found . 


We have used the recent reanalysis of p meson distribution amplitudes to improve 
and update LC sum rules for the semileptonic form factors, including hrst estimates of 
higher-twist corrections. Our hnal results for the form factors, decay rates and the spectra 
are presented in Tables 1,2 and in Figs. 7,8 together with lattice data and the results of 
quark models. We have given a detailed analysis of uncertainties of our approach, with the 
conclusion that its present accuracy is sufficient for a model-independent determination of 
\Vub\ with an error less than 20%. 

The accuracy of our results can be improved, by calculating radiative corrections to the 
sum rules and higher-twist corrections. Both is possible using existing methods and could 
ultimately decrease the uncertainty by a factor two, of order ~ 10% in |14fe|- Yet higher 
accuracy is, however, not feasible within the sum rule method. 

Acknowledgements: P.B. is grateful to the theory group of NORDITA in Copenhagen 
for its hospitality while this work was hnalized. 


Note added: When this paper was in writing, the work appeared with a LC sum rule 
for In the SU(3) limit their formula agrees with ours (except for the ^-function 

terms related to continuum subtraction). 


Appendices 


A Continuum Subtraction in LC Sum Rules 


The “standard” procedure, to which we conform in this paper, consists in approximating 
the (unknown) physical spectral function by the perturbative one above some threshold 
So, so that 



pPhyB(5 

s — 



^pert(g) 

S — 


Thus it is necessary to know the perturbative spectral function explicitly. 

In evaluating the correlation function (|2.12|) one encounters terms of type (g = ps —Pp, 
f{u) arbitrary function): 


In 



duf{u) 



1 

{q + uppYY' 


(A.l) 










The dispersive representation of Ji is trivial and reads 

poc /*! 


Ji = — / ds ■ 


/ du f{u) 6{us — ml — uuml + tu). (A.2) 

Jo ^ 


{Pp + qy 

Putting the upper limit of integration in s to Sq simply introduces a factor 0[c(m, s^)] 
in the integration over u. The function c is dehned after Eq. ( p.26|) . For higher powers 
one has either to integrate over u by parts, or calculate the spectral function by applying 
consequently two Borel transformations p\ —>■ and 1/M^ —> s, see [^, for details. 

In particular, for I 2 we hnd 

^ ^ 


Phis) = 


[mt 


t + 




mt 


t + ulm^^ 


(A.3) 


where m* is the solution of c(m*,s) = 0 inside the interval 0 < m* < 1. With this spectral 
density, performing the continuum subtraction and the Borel transformation in one 
obtains after a suitable change of variables 


Jo — continuum = 


M2 




/no 




uM^ 


{ml — ut + uum?) 


fin) 


mi — t + umt 


exp 


mM2 


[m. 


ut + uum^) 


Uq 


(A.4) 


where Uq is the solution of c(mo! Sq) = 0 with 0 < Mq < 1- Since in our case f{u) vanishes 
at M = 1, one arrives at the typical structure Q/{uM^) + 6 that enters (|2.25|) and ( p.26|) . 


B 3pt Sum Rules with Tensor-Current 

In this appendix we give the Wilson-coefficients entering the OPE of Eq. (|3.21|) . We 
use the invariant decomposition 

- g'^X) + ^+iPB + Pp)^Wp - «) 


+ {e^'apPsPpPp - (^^\pPBPpPi) +■■■, (B.l) 

where {Eq , E+, Sy } determines the form factors {Ai, A 2 , E}. Taking into account pertur¬ 
bation theory, the quark and the mixed condensate, as well as the four quark condensate 
(in vacuum saturation approximation), the OPE reads: 

s = (qq) + (qgaGq) + a,7r(qqy + ... (B.2) 

We give explicit formulas for the Borelized expressions and the double spectral 

function of Sp®’’*, such that 

yjpert _ f ^ -1- subtractions. 

J iSb-PB)iSu-pl) 
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87r2AV2 + 87r2A3/2 2sfeS„)(T 26), 


3b 


pT = 


{bT - 2sbSu) 


4^2A3/2 4^2^572 

+ 2sbSu{b‘^ + 2bsu + 3sf,Stj)}, 


{ 6(6 + 2su)T‘^ — 3su{b‘^ + 2bsb + SbSu)T 


pr = 


(6T - 2sfeS„) - {6^T^ - 66 s 6S„T + 2sbS„(6^ + 3sfeS„)}, (B.3) 


27r2A3/2 


27r2A5/2 


with A = sf + s2 + 6 — 2s6St, — 2sbt — 2s„6, b = Sb — ml and T = Sb + Su — t. For the 
nonperturbative terms we obtain: 


A3) 


s^W) _ 

Zjn — 


MlM^p 




sf = 0 , 
s® = 0 , 


^(5) _ mb mliMl ( rrij-t ^_ 
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3M^ 3M2 18M^ 36M| ISM^M^ 
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